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1. Introduction 

The magnetic moment and magnetic polarisability are fundamental properties of a particle 
which describe its response to the application of an external magnetic field. These properties can 
be calculated from first principles using lattice QCD techniques. One method is to use three-point 
functions with an explicit current in order to calculate electromagnetic form factors, from which 
the moment can be derived [1, 2]. Here we present results from the background field method [3- 
7] , which involves adding a phase factor to the gauge field which induces a constant magnetic field 
across the whole lattice. This field causes a shift in the energy of the particle which can be measured 
and used to calculate both the magnetic moment and magnetic polarisability via the energy-field 
relation 

E{B)=MN + li-B+^^-^^B^ + i^{B^). (1.1) 

The phase factor used to apply the magnetic field has the exponential form &xp{iaqA^{x)). 
This leads to a quantisation condition on the available choices for the magnitude of the magnetic 
field. If the size of the lattice is too small then the field will be large and the energy-field relation 
will no longer be dominated by the leading terms. Previous calculations using this technique have 
avoided the problem of the quantisation condition by using a linearised version of the phase factor 
with Dirichlet boundary conditions [5, 6], however this introduces new finite volume errors. The 
results presented here are calculated with the full exponential phase factor on a periodic lattice. 
These results include the magnetic moment and magnetic polarisability of the neutron, as well as a 
prehminary look at the magnetic moment of the lowest lying negative parity nucleon states. 



2. The Background Field Method 

To add a background magnetic field to lattice calculations we begin by considering the con- 
tinuum case, modifying the covariant derivative with the addition of a minimal electromagnetic 
coupling, 

D^=d^,+gGf,+qA^, (2.1) 

where is the electromagnetic four-potential and q is the charge on the fermion field. On the 
lattice this is equivalent to multiplying the usual gauge links by a simple phase factor 

u'j^\x)=QMiaqA^{x)). (2.2) 

To obtain a uniform magnetic field along the z-axis we note that B = V x A, and hence 

B, = dAy-dyA^- (2.3) 

Note that this equation does not specify the gauge potential uniquely. There are multiple valid 
choices of A^ that give rise to the same field, however they are equivalent up to a gauge trans- 
formation. We choose Ax{x,y) = By to produce a constant magnetic field of magnitude B in the z 
direction. Examining a single plaquette in the (/x, v) = {x,y) plane shows that this gives the desired 
field for a general point on the lattice. However on a finite lattice (0 < x < A^^v — 1 ) i (0 < j < A',. — 1 ) 
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there is a discontinuity at the boundary due to the periodic boundary conditions. In order to fix this 
problem we make use of the d^Ay term from Eq. (2.3), giving Ay the following values, 

Ay{x,y) = { ' (2.4) 

y-NyBx foiy=Ny-\, 

such that we now get the required value at the y =Ny — \ boundary. 

This still leaves the double boundary, x = Nx — l and y = Ny—\, where the plaquette only has 
the required value under the condition exp{—ia^qBNxNy) = 1. This gives rise to the quantisation 
condition which limits the choices of magnetic field strength based on the lattice size 

9 271?! 

where n is an integer specifying the field strength in multiples of the minimum field strength quan- 
tum. 



3. Simulation Details 

The results presented here were calculated on 2+1 flavour dynamical-QCD gauge field config- 
urations on 32^ X 40 lattices. These were provided by the PACS-CS collaboration [8] as part of the 
International Lattice Data Grid (ILDG) [9]. They have a physical lattice spacing of a = 0.0907 fm 
and a beta value of /3 = 1 .9. There ai^e five values of the light quark mass corresponding to the pion 
masses m^^ = 622, 512, 388, 282, 151 MeV, however the lightest mass has not yet been consid- 
ered herein. It should be noted that the configurations are dynamical only in the QCD sense, they 
have not been calculated with the background field included as this is prohibitively expensive. The 
corrections due to this are expected to be small. Unless mentioned otherwise correlation functions 
were calculated using the interpolating field Xi = {u^CYsd)u with 100 sweeps of stout smearing at 
the source. 



4. Magnetic Moment 

Once the background field method has been used to calculate correlation functions with a 
magnetic field on the lattice we use the energy relation of Eq. (1.1) to determine the magnetic 
moment. The term we want to isolate is /I • B, which we can do by taking advantage of the fact that 
the sign of this term is changed when the spin of the particle is flipped, while the magnitude stays 
the same. We simply take the difference of spin-up and spin-down energies, 

1 

r 

which is done by constructing the following ratio of correlation functions. 



-XE^-E{) = ix-B, (4.1) 



leaving us with only the term we are interested in. Although including the zero field correlators is 
not strictly necessary, it is useful for allowing correlated en^ors to cancel as much as possible. 
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Neutron Magnetic Moment 
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Figure 1: Left: Plot of the field strength vs the spin-difference mass shift. The solid line is a purely linear 
fit to the first two point and the dashed line is a linear plus cubic fit to all four points. Right: Plot of the 
magnetic moment vs pion mass squared showing a comparison of these results with older, quenched results. 
The left-most point marks the experimental value [10]. 



Figure 1 shows fits of the spin-difference mass shift to the field strength. These fits give the 
value of the magnetic moment from the linear coefficient. Both a purely linear fit to only the first 
two field strengths and a linear plus cubic term fit to all field strengths were calculated. This was 
to ensure that the fields used are small enough for the energy field relation (1.1) to be valid. Since 
the linear terms from each fit agree well within errors we can be confident that higher order terms 
are not unduly affecting the result. 

Figure 1 also gives the plot of the neutron magnetic moment results. Comparing quenched and 
dynamical values we see that they agree within eiTors and are both trending generally towards the 
physical value. Once the lightest mass has been calculated the trend should become clear and it 
will be relatively easy to do a chiral extrapolation to the physical pion mass and infinite volume. 

4.1 Negative Parity Nucleon Magnetic Moment 

It is also possible to examine the lowest lying negative parity nucleon state using the same 
correlation functions as for positive parity. This can be seen by writing the correlation function as 
a sum over positive and negative parity baryon states, 

(~'V,p) = l,A.B+AB+e yAs-AB-e s . (4.3) 

B+ '^'^B+ B- '^'^B- 

Then after setting 7? = to give Eb± = Mb± the negative parity spin-up and spin-down parts can 
be projected into the (3,3) and (4,4) components of the Dirac matrix via the parity projector, Y± = 
(l±7o)/2. 

Effective mass plots for the standard X\ type interpolating field showed poor signal. There- 
fore we decided to try using the interpolating field type Xi = {^i^ Cq)ysq. Figure 2 shows spin- 
difference effective mass plots for the neutron and the proton using a X2 interpolating field. Each 
one shows clear signal and an obvious linear progression over the fields as expected from the mag- 
netic moment effect. 

Whilst these plots appear to give fairly clear results for the magnetic moment, we must be 
careful. The two lowest lying negative parity nucleon states are the the N(1535) and the N(1650). 
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Figure 2: Spin-difference effective mass plots for the negative parity proton and neutron using Xi interpo- 
lating fields at the heaviest quark mass with all four fields strengths. 



The difference in mass between them is quite small, therefore we expect a lot of mixing. This means 
we cannot designate a specific state to each of the effective mass plots. In order to differentiate these 
states a variational analysis of a correlation matrix will be required, after which a real result for the 
negative parity nucleon magnetic moments can be reported. 

5. Magnetic Polarisability 

Similar to the magnetic moment case, calculating the magnetic polarisability requires isolating 
the desired term in the energy relation of Eq. (1.1). This time we take the average of spin-up and 
spin-down and also subtract the zero field energy in order to remove the bare mass, 

i((£t-Mt) + (£;-M;)) = ^-^/3B2, (5.1) 
which in terms of correlation functions looks like. 

In addition to the required magnetic polarisability term there is also another term left over. This 
term is the Landau energy and cannot be isolated from the polarisability term. Fortunately this term 
is zero for the neutron due to its lack of overall charge, so it can be ignored in this case. 

Even without any Landau energy to worry about the magnetic polarisability is significantly 
harder to extract than the moment. First of all the effect is much smaller at the important field 
strengths. Secondly the energy shift due to the polarisability arises because the presence of the 
field deforms the hadron away from its usual spherical shape. Since it is typical to use spherically 
smeared sources the creation operators have poor overlap with polarised states and it takes time for 
the deformation to appear. As seen in the left hand plot of Fig. 3 it is possible for the signal to start 
to become noisy before the energy reaches a proper plateau, making it difficult to fit. 

In order to get better plateau behaviour we examined different levels of smearing, looking at 
100, 35 and 16 sweeps of smearing. We also tried a point source in the hope that removing all bias 
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Figure 3: Left: The spin-averaged effective mass shifts for the heaviest quark mass at all four field strengths. 
Right: Quadratic plus quartic fits of the mass shift to the field strength. 
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Figure 4: The shift in the effective mass due to the field for spin-up and spin-down at the heaviest quark 
mass and smallest non-zero field strength. For spin-up the lines are, from top to bottom, 100 sweeps of 
smearing, 35 sweeps, 16 sweeps and zero sweeps. For spin-down the order is reversed. 



towards any shape would allow it to conform more easily to the final shape of a polarised neutron. 
Figure 4 shows plots of the energy shift due to the field for spin-up and spin-down. For each spin 
the point source has significant excited state contamination, leading to poor plateaus. The smeared 
sources all give similar results at large time, indicating the isolation of the ground state just prior 
to signal loss. 

It can be seen that the best plateau behaviour for spin-up comes from the 16 sweeps of smear- 
ing, whilst the best behaviour for spin-down comes from 100 sweeps. In order to take advantage 
of this difference we combined the spin-up 16 sweeps of smearing correlation functions with the 
spin-down 100 sweeps of smearing ones to create a new spin-average. Whilst the combination had 
slightly better plateau behaviour the improvement was not very significant, but it may be possible 
to improve the result further using a variational analysis of multiple smearings in order to pick out 
the best parts of each one. Another approach would be to use anisotropic smearing to try and match 
the shape of the polarised wavefunction, however this would likely require significant tuning. 
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6. Conclusion 

We have presented a variety of calculations of magnetic properties of the neutron using the 
background field method with a uniform background field. Results for the magnetic moment are 
clear and precise. There is good agreement between these dynamical results and previous quenched 
results and a good approach to the physical value. Magnetic polarisability results are more difficult 
but still promising. We have examined ways in which to improve the plateaus needed to fit for 
the polarisability and hope to improve these further using more advanced techniques. We have 
performed preliminary calculations of the magnetic moment of negative parity nucleon states. We 
found strong signal using the X2 interpolating field but require a more detailed analysis in order 
separate the two lowest lying negative parity states. Overall the background field method continues 
to be a valuable tool for investigating the magnetic properties of hadrons. 
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